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The sensitivity of quantum magnetometers [1] is challenged by control errors and, especially in the 
solid-state, by their short coherence times. Refocusing techniques can overcome these limitations 
and improve the sensitivity to periodic fields, but they come at the cost of reduced bandwidth and 
cannot be applied to sense static (DC) or aperiodic fields. Here we experimentally demonstrate 
that continuous driving of the sensor spin by a composite pulse [2] known as rotary-echo (RE) [3] 
yields a flexible magnetometry scheme, mitigating both driving power imperfections and decoher- 
ence. A suitable choice of RE parameters compensates for different scenarios of noise strength 
and origin. The method can be applied to nanoscale sensing in variable environments or to realize 
noise spectroscopy. In a room-temperature implementation based on a single electronic spin in dia- 
mond [4, 5], composite-pulse magnetometry provides a tunable trade-off between sensitivities in the 
/iTHz~' range, comparable to those obtained with Ramsey spectroscopy [5], and coherence times 
approaching T\ [6]. 



Solid-state quantum sensors attract much attention 
given their potential for high sensitivity and nano appli- 
cations. In particular, the electronic spin of the nitrogen- 
vacancy (NV) color center in diamond is a robust quan- 
tum sensor due to a combination of highly desirable prop- 
erties: optical initialization and read-out, long coherence 
times at room temperature (Ti > 2 ms [7, 8], T 2 > 0.5ms 
[9]), the potential to harness the surrounding spin bath 
for memory and sensitivity enhancement [10], and bio- 
compatibility [11]. 

Magnetometry schemes based on quantum spin probes 
(qubits) usually measure the detuning Suj from a known 
resonance. The most widely used method is Ramsey 
spectroscopy [12], which measures the relative phase Suot 
the qubit acquires when evolving freely after prepara- 
tion in a superposition state. In the solid state, a severe 
drawback of this scheme is the short free-evolution de- 
phasing time, Tg, which limits the interrogation time. 
Dynamical decoupling (DD) techniques such as Halm- 
echo [13] or CPMG [14] sequences can extend the coher- 
ence time. Unfortunately, such schemes also refocus the 
effects of static magnetic fields and are thus not applica- 
ble for DC-magnetometry. Even if Suj oscillates with a 
known frequency (AC-magnetometry), DD schemes im- 
pose severe restrictions on the bandwidth, as the opti- 
mal sensitivity is reached only if the field period matches 
the DD cycle time [1]. Schemes based on continuous 
driving are thus of special interest for metrology in the 
solid-state because they can lead to extended coherence 
times [15]. Recently, DC-magnetometry based on Rabi 
frequency beats was demonstrated [6]; in that method, a 
small detuning along the static magnetic field produces 
a shift « of the bare Rabi frequency ft. Despite ide- 
ally allowing for interrogation times approaching T\ , lim- 
iting factors such as noise in the driving field [6] and the 
bad scaling in Suj <C £1 make Rabi-beat magnetometry 
unattractive. More complex driving modulations [2] can 



provide not only a better refocusing of driving field inho- 
mogeneities, but also different scalings with Suj, yielding 
improved magnetometry. In this work, we use a novel 
composite-pulse magnetometry method as a means of 
both extending coherence times as expected by continu- 
ous excitation, and keeping a good scaling with Suj, which 
increases sensitivity. The i?-RE is a simple composite 
pulse (Fig. l.a) designed to correct for inhomogeneities 
in the excitation field; here, $ parametrizes the rotation 
angle of the half-echo pulse. For ^ 27rfc, fceZ, RE docs 
not refocus magnetic fields along the qubit quantization 
axes and can therefore be used for DC-magnetometry. 
For i? = 27rfc, RE provides superior decoupling from both 
dephasing [16] and microwave noise and can be used to 
achieve AC-magnetometry. 



In the rotating frame associated with the microwave 
field, and applying the rotating wave approximation, the 
Hamiltonian describing a continuous stream of i9-REs is 

U(t) = - [CI SW(t)a x + Su(l- a z )] , 

where 8W(t) = ±1 is the square wave of period T = 



On resonance (Suj = 0) the evolution is governed by the 
propagator U = e i T TW (*) CT - ; with TW(t) the triangular 
wave representing the integral of SW(t). We approxi- 
mate the time evolution in the presence of a detuning Suj 
by a first order Average Hamiltonian expansion [17, 18], 
yielding an effective Hamiltonian over the cycle, 
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Extending the above approximation to include the fast 
Rabi- like oscillations of frequency ($^\ 2 n) > wc can thus 
calculate the population evolution for one of the qubit 
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states, 
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The signal S indicates the presence of two spectral lines 
at M ± 2J f «n (!) for a detuning 5u. 

Thanks to this linear dependence on Suj, we expect a 
favorable scaling of the sensitivity rj, given by the shot- 
noise-limited magnetic field resolution per unit measure- 
ment time [1, 19]. For N measurements and a signal 
standard deviation A5, the sensitivity is 
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where 7 e (w 2.8MHzG^ 1 for NV) is the sensor gyromag- 
netic ratio and AB is the minimum detectable field. We 
broke down the total measurement time T into interro- 
gation time t and dead-time td required for initialization 
and readout. In the absence of relaxation, and neglect- 
ing td, a RE magnetometer interrogated at complete echo 
cycles t = n^f yields mE = ^ 2sin ^/ 2 ) - As shown 
in Fig. l.b, RE magnetometry has thus sensitivities com- 
parable to Ramsey spectroscopy, rm am sa ~~ ^f" Con- 
versely, Rabi-bcat magnetometry has r\R a bi ~ a t 
large times [18], which makes it unsuitable for magne- 
tometry despite long coherence times. 

To establish the sensitivity limits of RE magnetometry 
and compare them to other DC-magnetometry strategies, 
we carried out proof-of-principlc experiments in single 
NV centers in a bulk electronic-grade diamond sample. 
A static magnetic field B« pa 100 G effectively singles out 
a qubit {|0) , |1)} from the NV ground-state spin triplet, 
as the Zecman shift lifts the degeneracy between the |±1) 
levels. The qubit is coupled to the spin-1 14 N nucleus 
that composes the defect by an isotropic hyperfine in- 
teraction of strength A pa 2ir x 2.17MHz. After optical 
polarization into state |0), we apply a stream of n RE 
cycles using microwaves with frequency ui close to the 
qubit resonance wo = A + j e B\\, where A = 2.87GHz is 
the NV zero-field splitting. Because of the hyperfine cou- 
pling, wo is the resonant frequency only when the nuclear 
state is m; = 0. At room temperature, the nitrogen nu- 
cleus is unpolarizcd and, while its state does not change 
over one experimental run, in the course of the N ~ 10 6 
experimental realizations, pa 2/3 of the times the qubit 
is off- resonantly driven by \5cj\ = A. A typical n-cycle 
RE fluorescence signal is plotted in Fig. l.c for $ = it 
and ft pa 27rx 17MHz, while to determine the frequency 
content of the signal we plot the periodogram in Fig. 2. 

The number of distinguishable frequencies increases 
with interrogation time at the expense of signal-to-noise 



ratio. RE-magnctomctry not only discriminates the fre- 
quency shifts due to A pa 2tt x (2.14 ± 0.03)MHz but, for 
interrogation times as short as 5 [is, also those due to a 
small residual detuning b pa 2ir x (0.17 ± 0.02)MHz from 
the presumed resonance. In contrast, under the same 
experimental conditions Rabi magnetometry does not 
discern such a detuning before an interrogation time 
pa 188^s [18]. With longer interrogation times ~ 15/xs 
as in Fig. 2.b, we can detect a frequency as small as 
6 pa 2ttx(64± 12) kHz. 

To determine the experimental sensitivities, we esti- 
mate by driving the qubit with varying w, at fixed 
interrogation times t (Fig. 3. a). For each t, in Fig. 3.b we 
plot the minimum — t-^tt VNt and compare it to the ad- 

justed theoretical sensitivity tj/(CxCa)- Here (CxCa), 
pa (5.9 ± 1.4) x 10~ 3 in our setup, is a factor taking into 
account readout inefficiencies and a correction for the 
presence of the hyperfine interaction [1, 18]. The sen- 
sitivities agree with the theoretical model, with optimal 
~ 10/iTHz~5. Improved sensitivities arc expected from 
isotopically pure diamond [20]; adequate choice of inter- 
rogation times or polarization of the nuclear spin can 
easily set Ca = 1, while C can be improved by efficient 
photon collection [21] or using repeated readout meth- 
ods [18, 22]. 

The sensitivity of a NV magnetometer is ultimately 
limited by the interaction of the quantum probe with the 
nuclear spin bath. We model the effect of the spin bath 
by a classical noise source along a z [23], described by 
an Ornstein-Ulhenbeck (OU) process of strength a and 
correlation time r c . In the limit of long r c (static bath), 
the dephasing time associated with RE (Ramsey) mag- 
netometry is T RE = g ^^ /2)| (T 2 * = T Ram = 
respectively (see [18] for the general case). While at 
the optimum interrogation time -y one has tire / "HRam = 
\/ 2sin(tf/ ^ - > magnetometry allows a greater flex- 

ibility in choosing the effective coherence time, as larger $ 
increase the resilience to bath noise. Thus one can match 
the RE interrogation time to the duration of the field one 
wants to measure. In addition, RE can yield an overall 
advantage when taking into consideration the dead-time 
td- If td S> T' Ram , as in repeated readout methods [22], a 
gain in sensitivity can be reached by exploiting the longer 
interrogation times enabled by RE magnetometry [18]. 
An even larger advantage is given by AC-magnetomctry 
with 27rfc-RE [24], since RE provides better protection 
than pulsed DD schemes [16, 24]. 

Excitation field instabilities along a x also accelerate 
the decay of RE and Rabi signals. However, provided 
the echo period is shorter than t c , RE magnetometry 
corrects for stochastic noise in Rabi frequency [18]. This 
protection was demonstrated experimentally by apply- 
ing static and OU noise (r c sa 200ns) in the excitation 
microwave, both with strength 0.0517. The results for 
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Rabi and {ir, 57r}-RE sequences in Fig. 4 clearly show 
that whereas the Rabi signal decays within 0.5/xs, 5tt- 
RE refocuses static excitation noise and presents only a 
very weak decay under finite-correlation noise after much 
longer interrogation times 3/x, in agreement with the- 
ory detailed in [18]; 7T-RE is robust against the same noise 
profiles. 

The unique ability of the RE magnetometer to adjust 
its response to distinct noise sources is relevant when the 
sample producing the magnetic field of interest is im- 
mersed in a realistic environment; moreover, the field 
source might itself have a finite duration or duty cy- 
cle. The advantage is two-fold: the protection from noise 
can be tuned by changing the echo angle, thus allow- 
ing the interrogation times to be varied. Techniques for 
repeated readout in the presence of a strong magnetic 
field > 1000G [18, 22] can at once improve sensitivities 
and enable the use of much lower qubit resonance fre- 
quencies ^MHz, preferable in biological settings. Addi- 
tionally, a RE-magnetometer can discriminate magnetic 
noise sources given the sensor's well-understood dccoher- 
ence behavior under different noise profiles, effectively en- 
abling noise spectroscopy for both <j z and er^-type noises. 

NV center-based RE magnctomctry could find use- 
ful application, for example, to sense the activity of 
differently-sized calcium signaling domains in living cells, 
more specifically in neurons. Transient calcium fluxes 
regulate a myriad of cell reactions [42]. The signaling 
specificity of such fluxes is determined by their duration 
and mean travelled distance between membrane channel 
and cytoplasm receptor. The smaller, faster-signaling 
domains have resisted thorough investigation via both 
diffraction-limited optical microscopy [42], and the use 
of fluorescing dyes, which do not respond fast or accu- 
rately enough to Ca 2+ transients [26] . The magnetic field 
produced by as few as 10 5 Ca 2+ , being diffused within 
~ 10/iS through a ~ 200nm domain, can be picked up by 
a nanodiamond scanning sensor [27, 28] with sensitivity 
~ 10^THz~2 placed at close proximity (~ lOnm) [18]. 
The trade-off between sensitivity and optimal interroga- 
tion time under RE magnetometry can be optimized to 
the characteristics of the signaling domain under study 
by a suitable choice of i?. 

In conclusion, we have demonstrated a quantum mag- 
netometry scheme based on composite-pulses. Its key 
interest stems both from the continuous-excitation char- 
acter, offering superior performance for solid state sen- 
sors such as the NV center, and from the possibility of 
tuning the sensor's coherence time and sensitivity in the 
presence of variable or unknown sensing environments, to 
protect from or map noise sources. Current technology 
enables immediate implementation of such scheme at the 
nanoscalc. 



METHODS 

The NV center is a naturally occurring point defect in 
diamond, composed of a vacancy adjacent to a substitu- 
tional nitrogen in the carbon lattice. The ground state 
of the negatively charged NV center is a spin-triplet with 
zero-field splitting A = 2.87GHz between the m s = 
and m s — ±1 sub- levels. Coherent optical excitation at 
532nm promotes the quantum state of the defect non- 
resonantly to the first orbital excited state. While the 
m s = state mostly relaxes with phonon-mcdiated flu- 
orescent emission (ps 650 — 800nm), the m s = ±1 states 
have in addition an alternative, non-radiative decay mode 
to the m s = state via metastable singlet states. Due 
to this property, each ground state is distinguishable by 
monitoring the intensity of emitted photons during a 
short pulse of optical excitation. Additionally, contin- 
uous optical excitation polarizes the NV into the m s = 
state. We apply a magnetic field (w 100G) along a crys- 
tal axis (111) to lift the degeneracy between the m s = ±1 
states and drive an effective two- level system m s — {0, 1} 
at the resonant frequency [wq ~ 3.15GHz) obtained by 
continuous wave electron spin resonance and Ramsey 
fringe experiments. In each experimental run, we nor- 
malize the signal with respect to the reference counts 
from the m s = {0, 1} states, where the transfer to state 
m s = 1 is done by adiabatic passage. An arbitrary wave- 
form generator is used to shape the control microwave 
field, which is delivered by a striplinc in contact with the 
sample. A detailed setup description is found in [18]. 
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FIG. 1. RE magnetometry scheme and sensitivity, a) Experimental control sequence composed of n RE, in which 
the phase of the microwave field is switched by ir at every pulse of duration ^. b) Magnetometry sensitivities t]re of # = 
{3-7T/4, 7r, 57r}-RE sequences (green, blue, black), showing the tunability with the half-echo rotation angle. The sensitivity has 
its global minimum t]re ~ 1.38/\/i (comparable to Ramsey magnetometry, purple) for i9 « and consecutively increasing 
local minima for •& « (2k + l)n. A decrease in sensitivity is followed by an increase in coherence time, which can approach Ti 
as in Rabi-beat magnetometry (orange), whose sensitivity is limited by £2. Sensitivities are simulated in the presence of static 
bath noise using parameters from the fit depicted in c). c) A typical n = 55 cycles RE normalized fluorescence for $ = n and 
Q. ~ 2tvx 17MHz (blue); the modulation in the signal is due to the hyperfine interaction with the 14 N nucleus. The signal is 
filtered for even harmonics of 7rfi/($mod27r) [18] and then fitted to Eq. 1 modified to include decoherence induced by static 
bath noise (red). 
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FIG. 2. The periodogram identifies the frequency content of the signal, a) Periodogram for 7T-RE sequence for 
increasing interrogation times (thicker lines) . The periodogram is denned as the squared magnitude of the Fourier transform of 
the time signal. A pair of symmetric peaks about Q signals the existence of one detuning Sui. The number of resolved frequencies 
increases with time, at the expense of signal-to-noise ratio. After 5/^s of interrogation, we can estimate both the hyperfine 
interaction A ~ 2n x (2.14 ± 0.03)MHz and a small residual detuning from the presumed resonance, b~ 2tvx (0.17 ± 0.02)MHz. 
In this estimate, we correct for the real rotation angle § m 0.96-7T using the difference between the nominal and experimentally 
realized Rabi frequency (symmetry point in the spectrum). The uncertainty in the measurement is estimated taking into account 
the total interrogation time, the number of points in the time-domain signal, and the S/N, as detailed in [18]. Periodogram 
peaks can be tested for their statistical significance [18, 34]; we confirm that all 6 frequency peaks are considerably more 
significant than ap = 0.01 significance level (red), b) Innermost pair of frequency peaks arising from a b ~ 2nx (64 ± 12)kHz 
residual detuning in another experimental realization, for an interrogation time of 15/xs. 
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FIG. 3. Experimental sensitivity of RE magnetometry. a) RE signals at fixed interrogation times as a function of 
the detuning Suj from resonance, from which | J||j | is numerically calculated to obtain the sensitivity 77. With increasing 
interrogation times, the slopes initially increase, indicating an improvement in 77; the effect of decoherence for the longer 
interrogation times degrades the sensitivity, and the slopes smoothen accordingly. The different amplitude modulations are due 
to the three-frequencies in the signal, {b, A ± b}; polarizing the nuclear spin [30] would eliminate this modulation. From the 
fitted resonances for each curve (red diamonds), we estimate the true resonance to be at 0.09 ± 0.15MHz from the presumed 
resonance. Typical errors in the measurement are indicated (red errorbars). Interrogation times are chosen to coincide with 
minima of the sensitivity in the presence of the hyperfine interaction; in other words, Ca is at a local maximum at those 
times [18]. b) For each fixed interrogation time, the minimum sensitivity r\ within one oscillation period of the fitted oscillation 
frequency obtained in a), r = 2tsln j 1? / 2 ) ; j s plotted. The experimental points agree in trend with the theoretically expected 

sensitivities 77/ (C xCa) (solid curves), here corrected for the presence of static bath noise, r\ — > r\e { - t ^ TRE ^ > . T' RE was computed 
using a T 2 * ~ 2.19 ± 0.15/is fitting from a Ramsey decay experiment [18]. 
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FIG. 4. Experimental implementation of microwave frequency noise, against which RE is robust. 40 realizations 
of both static (red) and OU (r c = 200ns, blue) microwave noise of strength 0.05f2, f! ~ 2irx 19MHz, for a) Rabi, b) 5-7T-RE 
and c) 7T-RE sequences. We note that the correlation time of the driving field noise was purposely set much shorter than 
what usually observed in experiments in order to enhance the different behaviour of 7T-RE and 5-7T-RE. Typical errors in the 
measurement are indicated (black errorbars). a) We plot the peak of the Rabi fringes in the presence of static (red crosses) 
and OU (blue circles) microwave noise, which are in agreement with the expected theoretical decay detailed in [18] (solid blue 
line and solid red line, respectively). The oscillations at the tail of the signals are due to the finite number of experimental 
realizations. The peaks of a no-noise Rabi experiment are plotted for comparison (black diamonds) . b) The peak of the 5-7T-RE 
revivals are plotted in the presence of static (red crosses) and OU (blue circles) microwave noise; the peaks of the 57T-RE in 
the absence of microwave noise are also plotted for comparison (black diamonds). Whereas the echo virtually does not decay 
in the presence of static noise, under the effect of stochastic noise, the echo decays only weakly as stipulated by theory (solid 
blue line) described in [18]. c) The peak of the 7T-RE revivals are plotted in the presence of OU microwave noise (blue circles); 
virtually no decay is found if compared to the peaks of a no-noise experiment (black diamonds). 
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SETUP 

Experiments were run at room-temperature with single NV centers from an electronic grade single crystal plate 
([100] orientation, Element 6) with a substitutional nitrogen concentration < 5 ppb. The fluorescence of single 
NV centers is identified by a home-built confocal scanning microscope. The sample is mounted on a piezo stage 
(Nano-3D200, Mad City Labs). The excitation at 532nm is provided by a diode-pumped laser (Coherent Compass 
315M). and fluorescence in the phonon sideband (~ 650 — 800nm) is collected by a lOOx, NA = 1.3 oil immersion 
objective (Nikon Plan Fluor). The fluorescence photons are collected into a single-mode broadband fiber of NA = 0.12 
(Font Canada) and sent to a single-photon counting module (SPCM-AQRH-13-FC, Pcrkin Elmer) with acquisition 
time 100 or 200 ns. 

Laser pulses for polarization and detection are generated by an acousto-optic modulator with rise time < 7ns 
(1250C-848, Isomet). A signal generator (N5183A-520, Agilent) provides microwave fields to coherently manipulate 
the qubit. An arbitrary waveform generator at 1.2GS/s (AWG5014B, Tektronix) is employed to shape microwave 
pulses with the help of an I/Q mixer (IQ-0318L, Marki Microwave), and to time the whole experimental sequence. 
Microwaves are amplified (GT-1000A, Gigatronics) and subsequently delivered to the sample by a copper microstrip 
mounted on a printed circuit board, fabricated in MACOR to reduce losses. 

A static magnetic field is applied by a permanent magnet (BX0X0X0-N52, K&J Magnetics) mounted on a rotation 
stage, which at its turn is attached to a three-axis translation stage; this arrangement enables the adjustment of the 
magnetic field angle with respect to the sample. The magnetic field is aligned along a [111] axis by maximizing the 
Zeeman splitting in a CW ESR spectrum. 

DYNAMICS UNDER ROTARY ECHO SEQUENCE 

Consider a two-level system, |0) and |1), with resonance frequency ujq. The qubit is excited by radiation of 
frequency uj with associated Rabi frequency f2 and phase modulation <p(t), such that the magnetic field amplitude is 
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f2cos(wt + <p(t)). The Hamiltonian is then 



Hlab 







H cos(wi + <£>(£)) 



n cos(cj< + ip(t)) ujq 
In a frame rotating with the excitation field, the operator 



transforms the Hamiltonian to 



H 





ncos(Lut + (p(t))e v 



1 

e lun 



f2 cos(wi + <p(t))e 

LlJq — LU 



Applying the rotating wave approximation and setting set 5lj = loq — w, the Hamiltonian reads 

fe^W 
1 ge-^W 5 W 



(3) 



(4) 



(5) 



(6) 



One rotary echo (RE) is composed of two identical pulses of nominal rotation angle ?9 applied with excitation phases 
shifted by ir. Under a sequence of RE, e^*) = SW(£), with §W(£) the square wave of period T = equal to the RE 
cycle time: 



4 ^ 1 fkirnt\ 



sw ( t) = - ^ 

fe=l,odd v 7 



On resonance (Jw = 0) the evolution operator is trivially obtained 

cos(§TW(t)) -isin(|TW(t) 
-isin(§TW(£)) cos(§TW(£)) 



where TW(£) is the triangular wave representing the integral of SW(£) , 

fc=l,odd v 7 



( -TW(£)j l-isinf-TW(£) 



(7) 



(8) 



(9) 



Using Uq we make a transformation to the toggling frame of the microwave [31] to obtain the Hamiltonian W: 



u 



Slo ( 1 - cos(OTW(£)) isin(OTW(£)) \ _ Su 



2 V -«sin(OTW(£)) 1 + cos(OTW(£)) 



[1 - cos{QTW(t))a z - sin(f2TW(£))cr J/ 



(10) 



7-L is periodic with T and has a strength T8u <C 1, and can thus be analyzed with an average Hamiltonian expansion. 
In order to do so, we first express the elements of % in their Fourier series: 



eo s(ffl W(£)) = ^ + 2tfsintf £ ^^cos(^f); 

fc=l,odd V 7 
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sin(mrW(£)) = + 2d 

fc=l,odd 



(-l) fe ((-l) fe - cosi?) fknnt\ 
- cos 1 



d 2 - k 2 7T 2 



V J 



To first order then, 



-cos (0/2) £ sin (0/2) 
-i sin (0/2) cos (0/2) 
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For n rotary cycles, the propagator is approximated by Ure = e'^** 1 ' sa e mT ^ ( ' _ The population of a system initially 



prepared in 1 0) , under the action of U re , is described at full echo times by the signal 

' ■ (14) 



S(n) « \ 



I $\ 2 {#\ f^Sum (<d 
I - I cos I sin I - 



Extending the above approximation to include the fast Rabi-like oscillations of frequency ^^\ 2 tt) ' we obtain 



zA . , f$\ /2Sut /t?\\ / TT^t \ 



1 + cos 2 [ — I + sin 2 I ^- ) cos I — — sin I — I I cos 



2 J V 2 / V ^ V 2 / / V (i?mod27r) 







(15) 



indicating the presence of two spectral lines at ^J^tt) 2 T i snl (f ) ^ or eacrl existing detuning Sio. 

Our numerical studies suggest the existence of further signal components arising from higher frequency components 
in the Fourier expansion, which are not contemplated by the first-order approximation outlined above. Such compo- 
nents are oc cos sin (|)) cos f ^^j^y) an d oc cos f ^j^ja^j ) , fc e Z, thus being linked to split pairs of spectral 
lines around g , and single lines at ^g^y. 

PERIODOGRAM 

The periodogram is defined as the squared magnitude of the Fourier transform (FT) of the signal S(t) at times tj 
(j = 1, . . . , M), V = jj\ J2jLi Aj^ t] I 2 , where d, are the M data points [32]. 

Unlike the FT, the periodogram does provide bounds for frequency estimation from spectral analysis, besides 
being able to accommodate for noise profiles beyond static and white noise [32]. Take a simple sinusoidal signal 
dt = K cos(27r ft) + et, where e< is the added noise characterized by a (least informative) Gaussian probability 
distribution Normal(0,<7 2 ), with a in circular frequency units. To a is assigned Jeffrey's prior —, which indicates 
complete ignorance of this scale parameter. Under these conditions, the estimate frequency content of the signal is 
given by f est = f peak ± Sf, where f peak is the frequency of the periodogram peak, and 

6f = ^^=. (16) 

7T KtVM 

Here t is the total interrogation time for the M data points. Sf correctly takes into account the effect of both the 
interrogation duration t and the S/N = K ™ s = and is shown to correspond to the classical Cramer- Rao 

bound [33]. Note that Sf is in general smaller than the so-called Fourier limit, SfFi = The method is readily 
applicable to signals with multiple frequency content {/»}. 

In Supplementary Fig. 5, we compare typical experimental periodograms for 7T-RE, Ramsey and Rabi signals taken 
under the same conditions, for increasing interrogation times. The Ramsey periodogram, despite its lower signal 
intensity, clearly shows the 3 detunings {b + 2irx 5MHz, A± (b + 2ttx 5MHz)} present in the signal after l/xs; 7T-RE is 
sensitive to the residual detuning b ~ 2-7rx0.17MHz as explained in the main text after 5/is; finally, the Rabi sequence 
would only become sensitive to b after an interrogation time ~ 188/is, which is reflected in the broad single peak of 
the periodogram. 

We estimate the statistical significance level p of individual peaks [34] . Letting I m be the intensity of m-th largest 
ordinate among the total M in the periodogram, and calculating 

= ^ f v-^ m T ' (1^) 

the statistical significance of the m-th peak p m is approximated by 

p m w (M - (m - 1))(1 - T m ) M - m . (18) 

To determine Sf, we first estimate the S/N for each periodogram peak by dividing the peak area by the noise floor 
below the line of p = 0.01. 

The periodogram, if plotted over the full spectrum as in Supplementary Fig. 6, exhibits very high peaks corre- 
sponding to the even harmonics of which are present in the signal, but which are not taken into account by the 
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FIG. 5. Experimental periodograms for increasing interrogation times up to 5fis. a) Ramsey (5MHz detuned from 
the presumed resonance), b) 7T-RE and c) Rabi sequences. There is a trade-off between signal intensity and sensitivity to the 
detunings {b, A ± 6} in the signal. 
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FIG. 6. Periodogram over the full spectrum. Frequencies corresponding to the even harmonics of Q, which are present in 
the signal, but which are not contemplated by the first order of average Hamiltonian theory, can be clearly identified. In the 
inset, the signal peaks arising from the frequencies of interest are shown for interrogation times much longer than the dephasing 
time T' RE . 

first order of average Hamiltonian theory. In the inset, the peak structure originated from the detunings of interest is 
plotted for times much longer than the dephasing time T' RE . 



EXPERIMENTAL SENSITIVITY 



For a fixed interrogation time t, and scanning the detuning from resonance Su, we expect to observe the signal 



6{du)) oc cos I - sin I — ] ] = cos [oujt) 



(19) 



with t = f sin (|). 

In every experimental run, reference curves are acquired along with the signal S; they are noted TZq for the |0) state 
as obtained after laser polarization, and Hi for the |1) state as calibrated by adiabatic inversion. The signal is then 
normalized as 



S 



(20) 
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a) b) 




1 2 3 4 5 1 2 

Interrogation time (|is) Interrogation time (|xs) 



FIG. 7. Fitted parameters used to determine the experimental sensitivity, a) The fitted periods r (blue dots) are 
linear in t and agree well with the theoretically expected time 2t sln ^W 2 ) ( rec j line), b) The Ramsey signal (blue circles) is fitted 

to SRam = ki — [cos(<5wt) + cos(( J 4 + 5(jj)t) + cos((A — 5u)t)] e"^ 1 ^ 2 ' (red), with fitting parameters {fci, fe, Sui, A, TJ }. In 
particular, T 2 * ~ 2.19 ± 0.15/Lts. 



The standard deviation of the normalized signal is readily obtained 
AS 



(AK ) 



S-Tii 



(n - -Rif 

The sensitivity is calculated for the whole signal 

T)(5U)) 



S-Ki 



{Tlo-Tli) 2 Ho-ll! 



1 &S r— 

^VNt 



2 


1 


+ (AS) 2 


Ko-Kl 



7e 



dS 



(21) 



(22) 



for each fixed interrogation time t, we single out the minimum sensitivity rj(Su>) within one period of the fitted 
oscillation period r, depicted in Supplementary Fig. 7. a. 

The standard deviation for the sensitivity measurements is obtained by 



1 

An = — 

7e 



d5 



i 



ASVNt = — 

7e 



1 - 25 



1 



^/S(l-S) |£j 



ASVNt 



(23) 



Additionally, to every point in the plot there corresponds a factor C taking into account imperfect state detection [1, 
35]. While the theoretical signal S represents the population in the |0) state, measured from the observable M = |0)(0|, 
the experimental signal records photons emitted by both |0) and |1) states, so that the measurement operator is best 
experimentally described by M' = nrj|0}(0| +ni|l){l|. Here, {no,ni} are Poisson-distributed variables that indicate 
the number of collected photons; if perfect state discrimination were possible, n — > oo and n\ — s- 0. Including this 
effect, after n full echo cycles, the signal is modified to 



S'(n) « - 



(3no + ni + (no — n\) cosi?) + (no — rii — (no — ni) cost?) cos ( — — — sin ^— 



(24) 



We calculate the sensitivity for S' in the best-case scenario of minimum sensitivity given by the accumulated phase 
( 4S £ n sin (j)) = f , and note the existence of a factor C, with respect to the ideal sensitivity, t]m' = Vm/C: 



c- 




( — lino + 5ni) cosd 
2(n -ni) 2 + 2 



(n + ni) 
(n - ni) 2 



8n 



(n -ni) 2 sin 2 (^/2) 



(25) 



We use for no (ni) the mean photon number for the |0) (|1)) reference curve during each acquisition for different t. 
On average, nTf ~ 0.0022 ± 0.0003 and W - 0,0015 ± 0.0002. 

We also consider the fact that the signal S' has contributions from three detunings {Soj, A ± Sui}, where A is the 
hyperfine coupling between the NV center and spin-1 14 N nucleus; taking such detunings into account is, incidentally, 
fundamental for the choice of interrogation times: given the modulation imposed by the multiple frequencies in the 
signal, full echo times yielding a high signal amplitude arc preferred. A different strategy would be to polarize the 
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nuclear spin [36]. In order to compare the ideal sensitivity with the experimental one, in our experiments we need 
to introduce a further correction factor Ca, since the accumulated phase is only equal to the optimal ^ for the 
experimental realizations with mj = 0. We expect the sensitivity to become larger as ija = Vm'/Ca = Vm/(C x Ca), 
with 



C~A l = 



1 + 2 cos 



/ 2 At sin(i?/2) \ 
V * ) 



(26) 



In order to estimate C x Ca, we use the fitted value for A at each point (A ~ 2ir x (2.21 ± 0.07)MHz), the time t 
corresponding to the number of cycles at which the experimental point was taken, and a corrected $ ~ 0.984-7T that 
takes into account the real angle, given the experimental Rabi frequency, imposed by the duration of the echo half 
cycle, which can be controlled only up to the inverse of the AWG sample rate. A mean total correction factor of 
C x Ca ~ (5.9 ±1.4) x 10~ 3 is obtained for the set of points. The mean sensitivity curve (solid line) is expressed as 
the theoretically expected sensitivity in the absence of noise tjm, divided by C x Ca- Similarly, the lower (higher) 
bounds for the sensitivity are estimated by dividing the theoretical sensitivity by the maximum (minimum) C x Ca 
value in the set of points. 

The effect of decoherence is included in the plot using a fit for ~ 2.19 ± 0.15/xs from the Ramsey experiment 
shown in Supplementary Fig. 7.b. Assuming static Gaussian noise, we let t)a — > ^e^ T »'' 2 , with T' h 



T, i5 



RE — 2sin(i?/2) 



RABI-BEAT MAGNETOMETRY 



Rabi-beat magnetometry using a single solid-state qubit was recently demonstrated [6] . The scheme presupposes the 
existence of an absolute frequency standard against which one wishes to resolve a nearby frequency. For magnetometry 
purposes then, 

S = \{S Ra u(5u) - S Rabt (0)) , (27) 
where Su> denotes a detuning from the frequency standard. The sensitivity reads 



1 AS n _ V2^ / m _ 

??_ >5™o|^-| 7e y 2-2cos(tn)-tnsin(ifi) ' ^ ' 



?y is close to minima at t w (2k + yielding 



V2~n i I 2 



7. 



fen~— A/l + ^, (29) 



which tends to v2T^/7e for increasingly large interrogation times. 



EVOLUTION UNDER BATH NOISE 



In the presence of Gaussian static noise in the ^-direction with variance a 2 , the RE signal decays as 



(Sre) = 



1 



1 + cos 2 



( cos (^^sin e {t/T ' RE)2 



(30) 



where we define the dephasing time 



T> g 
RE <7>/2|sin(t?/2)|' 



(31) 
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Similarly, one obtains for the Ramsey signal 

1 



-(t/T' Ram y 



with T Rarn = T 2 = ~ 

a 



(32) 



Note that T' RE > T' Ram always; nevertheless, at the optimum interrogation time calculated for both sequences as 



Vre 

'} Ram 



2sin(iV2) 3 



> i ; 



(33) 



the sensitivity ratio above has a minimum ijre /VRam ~ 1-20 for i? ~ which is the angle that yields the highest 
sensitivity for the RE sequence. 

We now turn our attention to the evolution of the Rabi signal under Gaussian dcphasing noise. For Suj <C n, the 
Rabi signal is approximately 



S 



Rabi 



= l 



o- 



fl 2 + Suj 2 



i-i- 



"n 2 



sin 2 [-in 



Suj 2 

7T 



calculating the expected value (SR a bi) under the noise distribution yields 



(SRabi) — — 



cos(tn + arctan(to 2 /n)/2) a 2 ( cos(tfl + 3 arctan(to 2 /n)/2) \ 



(34) 



(35) 



In the presence of stochastic (Ornstein-Uhlenbeck) noise with zero mean and autocorrelation function a 2 e Tc , a 
Ramsey signal decays as [37] 



(S Ram ) = \ (l + e-C'W) , with C(t) = cr 2 T 2 c {t/r c + e~£ - 1) . 
Numerical simulations valid for t c <j < d/2 and r c > i)/(2fl) indicate that the RE signal decays as 



(Sre) = - 



( | 



with 



WW*) 



/J 2 



(36) 



(37) 



(38) 



Note the additional factor 



4sin 2 (i9/2) 
i3 2 



( CT ) 



< 1. 



Simulations that compare different i9-RE for $ = {37r/4, n, 5tt} and Ramsey signals in the presence of stochastic 
noise are depicted in Supplementary Fig. 8. a; it is clear that RE sequences are more resilient to bath noise with 
correlation times shorter than the echo period. 

Previous calculations [38] indicate that, for slow baths — -C the Rabi signal follows the static noise behaviour 

at 2 4fi 2 

for short times, and decays oc e 2 %^r for long times. Fast baths ^- 3> induce a decay of the Rabi signal cx c 4t <= . 



EVOLUTION UNDER EXCITATION FIELD NOISE 

In the presence of a constant error in the Rabi frequency such that S7 — > (l+e)n, the infidelity (1 — Tr[£/(e)C/(0)]/2) = 
(1 — F) of the pulse sequence is given to second order in the detuning from resonance Suj and in e by 

. e 2 t 2 8uj 2 (2 + tf 2 -2cos?9-2t?sintf) 

(1 - F) RE m (39) 
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FIG. 8. Signal decay in the presence of stochastic noise, a) In the presence of stochastic bath noise, different $-RE with 
"D = {37r/4, 7r, 5ir} (green, blue, black) are compared against a Ramsey sequence (purple); numerical simulations (dashed lines) 
agree with the formulas presented in the text (solid lines). The Ramsey sequence is the least resilient to bath noise, whereas 
one can adjust the dephasing of the RE by the choice of RE sequences are more resilient to bath noise with correlation times 
shorter than the echo period. The used numerical parameters are: Q, = 2ir x 20MHz, Sui = 27rx2MHz, t c = 200ns, a = 0.05Q,. 
b) In the presence of stochastic noise in the excitation field, the situation is inverted: RE sequences refocus microwave noise 
with correlation times longer than the echo period; the decay of the Rabi sequence (orange) is plotted for comparison. Used 
parameters are the same as in a), except for Slj — 0. 



for RE and 

e 2 t 2 n 2 e 2 Suj 2 (-2 + t 2 il 2 + 2cos(tfi)) , s 

(1 - F) RaU *— 8 * m ^ (40) 

for Rabi-beat magnetometry. 

Similarly, an error in the Rabi frequency will yield a flip-angle error in the Ramsey sequence, resulting in the 
infidelity 

e 2 7r 2 e 2 &j 2 (-16 + 4-7r 2 + 7rfO(8 + 7rffi)) 
(1 - F) Ram « — — 3 . (41) 

„ t_ 

In the presence of stochastic noise in the excitation field with zero mean and autocorrelation function a e T <= , the 
resonant cases for RE, Rabi have simple analytical solutions. 

A cumulant expansion technique applied to periodic Hamiltonians [37, 39] yields for the envelope of a resonant RE 
sequence 

(S RE ) = \(l + e-^) , (42) 



with 

"2m9 _ , _ . , /l ■d 



C(n) = r 2 c o 2 



(JTr 



+ 1) + e r c- 



(43) 



We note that this decay is equivalent to the decay under pure dephasing for a PDD sequence [40]. 

In Supplementary Fig. 8.b, we simulate the signal for different i9-RE and Rabi sequences if noise in the excitation 
field is present. Contrarily to RE decay in the presence of bath noise, and as shown experimentally in the main text, 
RE sequences can refocus excitation noise with correlation times longer than the echo period. 

We note that the Rabi signal decay for noise along a x should be comparable to Ramsey signal decay in the presence 
of stochastic noise along a z . We thus have the decay 

(Sn. abl ) = \ (l + , with C(t)=<j 2 T 2 (t/T c + e-^ -1) . (44) 



The advantage of the RE sequence over the Rabi is thus the same advantage that dynamical decoupling sequences 
can offer. 
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FIG. 9. Sensitivity with repeated readouts. We compare the achievable sensitivity for Ramsey (purple, dotted) and 
Rotary echo (-7T-RE, blue, dashed; 117T-RE, gray) sequences, when using the repeated readout scheme, with N r = 100 and the 
time per each readout t T = 1.5/is. In the presence of dephasing with TJ = 3/is, the longer angle RE achieves good sensitivity 
for much longer interrogation times. 



REPEATED READOUTS 



efficiency C of the NV center spin state is thus given by Eq. 25, which for d = kn reduces to C = ( 1 - 1 - 3 (" 0+ " 1 



The NV spin state can be read under non-resonant illumination at room temperature using the fact that the 
m s = ±1 excited states can decay into metastable states, which live for ^300 ns, while direct optical decay happens 
in about 12 ns. Thus, a NV in the m 8 = state will emit, and absorb, approximately 15 photons, compared to 
only a few for a NV in the m s = ±1 states, yielding state discrimination by fluorescence intensity. Unfortunately the 
metastable state decays primarily via spin-non conserving processes into the m s = state thereby re-orienting the 
spin. This is good for spin polarization, but erases the spin memory and reduces measurement contrast. The detection 

-1/2 

(no-ni) 2 J 

where m,o is the number of photons collected if the NV spin is in the m s = {0, 1} state, respectively. 

In the repeated readout scheme [22, 41], the state of the nuclear spin is repetitively mapped onto the electronic 
spin, which is then read out under laser illumination. The measurement projects the nuclear spin state into a mixed 
state, but the information about its population difference is preserved, under the assumption that the measurement 
is a good quantum non-demolition measurement. We can include the effect of these repeated readout by defining 

N r (no- 
number of measurements. The sensitivity needs of course to be further modified to take into account the increased 
measurement time. Provided the time needed for one measurement step is smaller than the interrogation time, 
it becomes advantageous to use repeated readouts. As RE increases the interrogation time, it can achieve better 
sensitivity than Ramsey magnetometry by using the repeated readout scheme (which is instead not advantageous for 
a simple Ramsey scheme), as depicted in Supplementary Fig. 9. 



a new detection efficiency, Cjv r = ( 1 + -rr- 3("o+"i) J . which shows an improvement oc \/N r , where N r is the 



CALCIUM SIGNALING DOMAINS 



Although virtually all neuronal reactions are regulated by diffusing Ca 2+ ions between membrane channel sources 
and cytoplasm target receptors, triggering specificity is ensured by the fact that such diffusion events are localized 
in time and space. The size of the signaling domain, understood as roughly the distance between channel and 
receptor (50nm to 0.5/mi), determines the diffusion timescale (/is to ms) and strength, the latter measured by Ca 2+ 
concentration (100 to 1/jM) [42, 43]. 

Let a flux with duration t and mean travelled distance d between membrane channel and cytoplasm receptor. The 
magnetic field at a distance r from a transient Ca +2 flux composed of I ions is estimated as 



with e the electron charge, and with the magnetic permeability of the cell approximated by /io, the vacuum perme- 
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ability. Therefore, the minimum required sensitivity to sense the afore-described calcium flux is 



r, [-= =V27r^-piV^. (46) 
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